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Mechanical detection of nuclear spin relaxation in a micron-size 
crystal 
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Abstract. A room temperature nuclear magnetic resonance force microscope (MRFM), fitted in a lTesla 
electromagnet, is used to measure the nuclear spin relaxation of 1 H in a micron-size (70ng) crystal of 
ammonium sulfate. NMR sequences, combining both pulsed and continuous wave r.f. fields, have allowed 
us to measure mechanically Ti and Ti, the transverse and longitudinal spin relaxation times. Because two 
spin species with different T\ values are measured in our 7/im thick crystal, magnetic resonance imaging 
of their spatial distribution inside the sample section are performed. To understand quantitatively the 
measured signal, we carefully study the influence of the spin-lattice relaxation and the non-adiabaticity of 
the c.w. sequence on the intensity and time dependence of the detected signal. 



PACS. 07.79.Pk Magnetic force microscopes - 76.60.-k 
87.61.Ff Magnetic resonance imaging: Instrumentation 



Nuclear magnetic resonance and relaxation 



1 Introduction 



For a long time research groups have looked for new ways 
of detecting electronic or nuclear paramagnetic resonance 
with better sensitivity. A review of different proposed meth- 
ods can be found in the introduction of Abragam's book 
[1]. In two seminal papers [2,3], Sidles recognized the ad- 
vantage of coupling the spin system to a mechanical oscil- 
lator for magnetic resonance imaging. In this technique, 
the force signal is proportional to the magnetic field gra- 
dient [4], which, in extremely inhomogeneous field, should 
allow high spatial resolution. The new technique is referred 
to as magnetic resonance force microscopy (MRFM) [5]. 

The first magnetic resonance force signal was detected 
by Rugar et al. in 1992 while exciting electron spin res- 
onance (eMRFM) in a 30ng crystal of diphenylpicrylhy- 
drazil [6]. Two years later, Rugar et al. reported the me- 
chanical detection of 1 H (protons) nuclear magnetic reso- 
nance (nMRFM) in 12ng of ammonium nitrate [7]. These 
two pioneering experiments demonstrate that a micro- 
fabricated cantilever, identical to the ones developed for 
atomic force microscopy, can detect the magnetic moment 
of a micron-size sample. In the case of nuclear magnetic 
resonance (NMR) [7], the achieved sensitivity of 10 13 spin, 



a Present address: Service de Physique de l'Etat Condense, 
CEA Orme des Merisiers, F-91191 Gif-Sur-Yvette, e-mail: 
oklein@cea.fr. 

b Present arfrfress.'Laboratoire de Physique des Solides, Uni- 
versite Paris-Sud, F-91405 Orsay. 



at room temperature and in a field of 2.4T, represents a 
substantial improvement over the standard coil detection. 

Significant progress was made in the past few years. 
In 1996, Zhang et al. mechanically detected the ferro- 
magnetic resonance (fMRFM) of yttrium iron garnet [8] . 
Imaging experiments with eMRFM [9,10], nMRFM [11, 
12] and fMRFM [13] were performed. A magnetic reso- 
nance torque signal in a homogeneous magnetic field [14] 
was also detected. Improvement of the force sensitivity by 
operating at low temperature [15-17] was demonstrated. 
Force maps of the sample were obtained with the magnetic 
probe placed on the mechanical resonator in eMRFM [18] 
and fMFRM [13]. The highest sensitivity reported so far is 
around 200 electron spin in a 1Hz bandwidth. The result 
was obtained by operating an eMRFM at 77K in a very 
large magnetic field gradient [19]. In 1996, Wago et al. 
demonstrated that a pulse sequence combined with fast 
adiabatic passages can allow to measure the nuclear spin- 
lattice relaxation time of 19 F in calcium fluoride at low 
temperatures [15]. The same method was used to mea- 
sure the longitudinal spin relaxation of H in ammonium 
sulfate at room temperature and normal pressure [20,21]. 
Recent eMRFM work in vitreous silica at 5K showed that 
the same principles can be also applied to study electron 
spin dynamics of E' centers with long T± [17]. 

In this paper, we report the first measurements of 
both the transverse and longitudinal nuclear spin dynam- 
ics of 1 H using mechanical detection. A very thin sam- 
ple is used to analyze if new phenomena might be spe- 
cific to small sizes. Our instrument is a simple home-built 
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Fig. 1. MRFMs are miniature version of the Faraday balance. 
The instrument fits between the poles of an electromagnet that 
generates a homogeneous static field B ext = IT. A sample is 
affixed onto a micro-fabricated cantilever and placed in a field 
gradient g — 500T/m produced by a 2 millimeter-size iron bar. 
The gradient serves to create a magnetic force on the spin sys- 
tem. The induced A-scale bending of the cantilever is measured 
by a laser beam deflection on a photodiode. To increase the 
sensitivity, the nuclear magnetization is inverted periodically 
at f.4kHz, a frequency that corresponds to the fundamental 
flexure mode of the cantilever. The induced vibration is moni- 
tored by a lock-in. The nuclear spin oscillations are obtained by 
cyclic adiabatic inversions using frequency modulated r.f. field 
produced by a VCO. The carrier is at the Larmor frequency of 
protons. 



MRFM located inside a lTesla electromagnet. The me- 
chanical motion of the cantilever is monitored by a laser 
beam deflection system. The sample is a 7fim thick crys- 
tal of (NrLt^SC^. Two spin-lattice relaxation times are 
observed, T\ s = 0.4s and Tu — 5s. The later value corre- 
sponds to the T\ reported in the literature for this com- 
pound [22,23]. The short relaxation, however, might be 
due to water contamination inside the crystal during its 
contact with air. These same two relaxation rates are also 
measured by conventional NMR in powder samples with 
particles of dimensions smaller than 50/im. 

After introducing in Section 2 the measurement tech- 
nique employed in this study, we will present in Section 
3 our results on the transverse and longitudinal spin re- 
laxation properties of (NH4)2S04. This will be followed 
in Section 4 by a more detailed analysis of the time de- 
pendence and magnitude of the force signal in order to 
quantify the properties of spin with short and long T\ and 
to determine the effect of the non-adiabaticity of the se- 
quence in the measured signal. Finally a model to describe 
our experimental data will be proposed. 



2 Measurement of the force signal 

The setup is schematically represented in Fig.l. The ex- 
periment is performed at room temperature inside a vac- 



uum cell (10 _2 torr) constantly connected to the inlet of a 
rotary pump. The instrument [24] fits between the poles of 
an iron core electromagnet which produces a static mag- 
netic field B ext k along the z axis. To the uniform field 
one adds a second inhomogeneous field with axial sym- 
metry produced by a magnetized iron bar 8mm in length 
and 1.9mm in diameter. The polarization field experienced 
by the spin is -Bo = B ext + B cy i, with Bq = Bq ■ k and 
B cy i — B cy i ■ k. Near the symmetry axis, the instanta- 
neous magnetic force acting on the sample is given by the 
expression [25]: 



F(t)= f M z (r,t) 
Jv, 



dB 



cyl 



dz 



dv. 



(i) 



Here M z is the z component of the bulk magnetization 
and V s is the volume of the sample. For small sample 
size, we make the approximation that the field gradient 
g = dB C yi/dz is uniform over V s . A new length variable 
( = Bo(r)/g is defined so that a plane of constant £ maps 
onto a surface (actually a paraboloid) of constant polar- 
ization field which also corresponds to a sheet where the 
spin have the same motion. A paraboloid of fixed C value, 
however, shifts axially away from the iron cylinder when 
B ext increases. In this experiment, the sample is placed 
0.70mm above the iron cylinder and centered on the cylin- 
der axis. At this distance, the calculated axial field gradi- 
ent is g = — 470T/m (see appendix A). 

The mechanical force detection is obtained by measur- 
ing the elastic deformation along the z axis of a micro- 
fabricated cantilever on which the sample is attached. In 
this orientation, the probe is sensitive to the longitudinal 
component of the nuclear magnetization in contrast with 
a standard coil detection. The cantilever equation of mo- 
tion is represented by a damped harmonic oscillator with 
a single degree of freedom. The measurement technique 
uses the optical deflection of a 4 /iW HeNe laser beam 
which reflects off the rear side of the cantilever onto a 
position-sensitive detector. 

Our test compound is (NtLj^SO,!. This non-magnetic 
insulator has a high proton density d = 6.4 x 10 22 1 H/cm 3 
and is in its paraelectric state above 223K. NMR mea- 
surements of the 1 H spin-lattice relaxation time at 300K 
in our powder [26] give T\ z ss 5s along the static field 
[27]. The 1 H linewidth is 5G and the second moment is 
Mi = 4G 2 at 295K [28]. Our sample is a crystal cleaved 
to a platelet aspect ratio and glued with epoxy on the 
end of a commercial SiaN4 amorphous cantilever of spring 
constant k =0.008N/m as can been seen in Fig. 2. Af- 
ter completing the assembly, the cantilever resonance fre- 
quency drops from 5.8kHz to 1.4kHz due to the sample 
mass [29]. The quality factor of the loaded cantilever is 
Q «4000 in vacuum. From the electron microscopy im- 
ages (Fig. 2), the sample dimensions are approximately 
100 x 50 x 7/im 3 with the smallest length (the thickness) 
oriented along the axial field. This represents a volume 
V a = 3.5 x 10~ 8 cm 3 or a mass m = 70ng and corresponds 
to N « 10 15 protons. The temperature of the cantilever 
holder is stabilized around +27°C during the measure- 
ment [30]. The nuclear magnetization at thermal equilib- 
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Fig. 2. Images of a commercial SisN4 (amorphous) cantilever: 
A 7/xm thick single-crystal of (NH^SCU sample is glued on 
the cantilever end with epoxy. The loaded cantilever has a res- 
onance frequency of 1.4kHz, a spring constant of 0.008N/m 
and a quality factor of 4000 in vacuum. The estimated sample 
volume is 3.5 x 10~ 8 cm 3 



rium is expressed by the Curie law Mq = (d^Bo/kBT) fc, 
with /i n = 1.4 x 10~ 26 J/T the proton magnetic moment, 
T = 300K and B = 1.3T the polarization field. This gives 
a magnetic moment M V S = 2.3 x lCr 16 J/T. 

In order to increase the sensitivity, M z is modulated 
at a frequency u> m close to u> c , the frequency of the fun- 
damental flexure mode of the cantilever. At the moment 
the optimal configuration uses cantilevers that have me- 
chanical resonance frequencies in the audio range and Lar- 
mor frequencies lvq which are several orders of magnitude 
larger (radio or microwave frequencies). Only two meth- 
ods have been used to create an oscillatory force on the 
cantilever: cyclic saturation and cyclic adiabatic inversion. 
They are restricted respectively to compounds that have 
spin-lattice relaxation times T\ either shorter or larger 
than the oscillation period of the cantilever. 

In our case, the modulation of M is generated by 
a continuous-wave (c.w.) sequence that consists of peri- 
odic adiabatic fast passages [1]. The radio- frequency (r.f.) 
source is a 35-75MHz Voltage Controlled Oscillator (VCO). 
The r.f. output field is amplified up to 7W and fed into 
an impedance matched resonating circuit (Q r f ~ 100) 
tuned to a fixed frequency, 54.7MHz. A small coil (3 turns, 
0.8mm in diameter) is in series with the tank circuit. 
The sample is 0.5mm away from this antenna. The nu- 
clear spin are irradiated for a few seconds by a linearly 
polarized r.f. field B x = 2Bi cos L u{t')dt' with ui(t) = 
i?sin(cj m t) + loq, a sine-wave modulation of the r.f. fre- 
quency around the proton Larmor frequency ujq = 7<?Co> 
where j/2ir = 4.258kHz/G is the nuclear gyromagnetic ra- 



tio. The surface of constant ( = ( is called the resonant 
sheet. The sinusoidal frequency modulation is started at 
a time t = 0. In a transformation to a rotating coordinate 
system with an instantaneous angular velocity u)(t)k, the 
apparent magnetic field is: 



B e (C,t)=B 1 i 



k. 



(2) 



is defined as the polar angle made by the apparent field 
with the external field. The magnetization, however, pro- 
cesses about the direction B e + 8/-fj, with 9 = d9/dt 
(see appendix B). A parameter for non-adiabaticity is de- 
fined with tana = 9/(j | B e |) the angle between the two 
vectors. Provided that the adiabatic condition a < 1 is 
satisfied, the spin system remains at all times in a state 
of internal equilibrium and M is parallel to B e as re- 
quired by Curie's law. The longitudinal magnetization is 
M z (C,t) =\ M | cos6>, where 



cos6> = 



gC-u(t)h 



^/{gC-^hf + Bf 



(3) 



For free spin, | M | is a constant of the motion [1]. This 
is no longer true in condensed matter because of spin- 
lattice relaxation. In our sample, however, the magnetiza- 
tion decay is slow compared to the modulation period. Un- 
der our measurement protocol, an extra defocusing orig- 
inates from the lack of adiabaticity of the modulation. 
In a first step, these effects are neglected and they will 
be considered in a more detailed analysis deferred to a 
later section (see also the appendix B). At time t = 0, 
Bi is assumed to be turned on adiabatically with the 
sample initially in thermal equilibrium. In this case the 
norm M reflects the state of the longitudinal magneti- 
zation immediately before the force measurement. Dur- 
ing the c.w. sequence, the oscillatory movement of M z {t) 
comes from the cos 6 factor. The value is expanded in time 
series cos# w ao + a\ sin(w TO i) [31] with ai the first har- 
monic Fourier component [11] (higher harmonics have a 
negligible effect on the motion of the cantilever). Because 
of the large field inhomogeneities, the amplitude of os- 
cillation depends on the location inside the sample. The 
resonant sheet, which is the paraboloid of constant Co, 
corresponds to the surface of maximum amplitude of os- 
cillation. The spatial dependence of ai(C) is the sensitivity 
profile, r is the half width at half maximum of this bell- 
shaped curve, r has the units of a distance and it defines 
the thickness of the slice probed. The amplitude of r de- 
pends on both fl and B\ [32]. The induced vibration is 
synchronously amplified by a lock-in technique through a 
single-pole low-pass linear filter of time constant r;. For 
uj m — lj c , the lock-in signal grows exponentially (an exact 
expression will be given in equation (11)) to the asymp- 
totic amplitude 



1 Qg 



V2 k X."»*(0«- 



(4) 



In conclusion, the maximum amplitude of vibration achie- 
ved by the cantilever is proportional to the longitudinal 
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Fig. 3. (a) Vibration amplitude of the cantilever measured by 
the lock-in for a (NH^SCU crystal containing 10 ° protons at 
300K in B ex t = 0.9425T. The trace corresponds to a single shot 
experiment with no averaging. The lock-in time constant is 
n = 0.3s. (b) Details of the start and end of the c.w. sequence. 
The crystal is irradiated for 3s by a r.f. field of Bi = 10G (up- 
per panel). The bottom panel shows the frequency waveform 
applied to the VCO around u)o/2n = 54.7MHz. The amplitude 
of the frequency modulation is f2/2n = 150kHz. The width of 
the curve is the digitalization noise of the oscilloscope (not the 
phase noise). 



magnetic moment inside the probed slice at the begin- 
ning of the c.w. sequence. In the ideal case of a uniform 
inversion of all spin inside the sample, the asymptotic am- 
plitude would be A tot = QgM V s /(kV2). 

In the upper panel of Fig. 3 the time dependence of 
the lock-in output A(t) is shown when the c.w. sequence 
is applied. The time delay between force measurements is 
set to 27s (> 5Ti) to ensure a steady state magnetization 
close to the thermal equilibrium value. The lock-in time 
constant is t; =0.3s which corresponds to an output noise 
of 4 A. The bottom panel of Fig. 3 displays the time de- 
pendence of B\ and oj at the beginning and end of the 
c.w. sequence. At the start, the amplitude of B\ is turned 
on from to 10G in 5ms when the frequency is well off- 
resonance, i.e. 400kHz below wo/27r. The frequency is then 
ramped to resonance in 7ms. Finally, the frequency mod- 
ulation of the r.f. field is applied for 3s with a deviation 
Q/2-k = 150kHz. For these settings, the calculated value 
of r = 7/xm is comparable to the sample thickness. 

Since the r.f. tank circuit is tuned to a fixed frequency 
the resonance is found by sweeping the external field B ext . 
There is no spurious vibrations of the cantilever induced 
by the r.f. fields when B ext is outside the resonance range. 
Fig. 3a shows the amplitude of the lock-in signal achieved 
in a one shot experiment at the resonance maximum, B ext — 



Fig. 4. The amplitude of the force signal (lock-in peak am- 
plitude averaged over Is around its maximum) is shown as a 
function of the width of a r.f. pulse applied 13ms before the c.w. 
sequence. Each point is the average of 16 c.w. sequences. The 
solid line is proportional to a damped cos <p with ip = "/t p Bi 
the nutation angle. The 6.4W r.f. power during the pulse cor- 
responds to a rotating field of Bi — 15G at the sample loca- 
tion. The settings for the c.w. sequence are a r.f. field of 10G, 
S7/2n = 50kHz and t; = 100ms. The inset is a schematic of 
the time dependence of Bi (solid line) and u> (dashed line) 



0.9425T. The maximum vibration amplitude is around 
40 A which corresponds to a signal to noise ratio of 20dB. 
The shape of the lock-in signal A{t) depends on the value 
B ext [17]. For 5 ext =0.9425T, i.e. Co set at the middle of 
the sample thickness, no steady-state vibrations of the 
cantilever are induced by the c.w. sequence and the lock-in 
signal decays toward zero for long sequence. On the other 
hand, for B ext ^ 0.9425T, an unbalanced partial repolar- 
ization of the magnetization occurs during each cycle and 
the lock-in signal decays to a finite value which changes 
sign for B ext smaller or larger than 0.9425T. 



3 Relaxation measurements 

In this section the nuclear spin dynamics of our sample is 
measured by applying a series of r.f. pulses before the c.w. 
sequence described above. 

In order to calibrate the strength of the r.f. field, a 
r.f. pulse of duration r p is applied, with an amplitude B\, 
13ms before the c.w. sequence. During this pulse, Mq ro- 
tates about B e through an angle <p. The angle obtained 
at the end of the pulse is ip = 7 | B e \ t p [15,17]. Within 
a few milli-seconds after the pulse, the nutated magneti- 
zation vector decays to its longitudinal component which 
then determines the amplitude of the maximum vibrations 
achieved by the cantilever during the force measurement. 
B\ is set at maximum power, i.e. 6.4W, for the pulse. The 
c.w. sequence uses a 2.9W r.f. field. Fig. 4 shows the lock- 
in output averaged over a Is time interval around its peak 
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Fig. 5. Measurement of the transient shape of the spin-echo: 
A -K/2-T a -iT pulse sequence is used to form a spin echo. The 
transverse magnetization is measured with the combination of 
a ir/2 pulse and the c.w. sequence. The amplitude of the force 
signal is shown as a function of r a + r b with a fixed r a — 17/us. 
A r.f. field of B\ = 15G is used for the pulses. The settings for 
the c.w. sequence are a r.f. field of 10G, Q/2n — 50kHz and 
n = 100ms. The solid line is the expected shape of the spin 
echo for a platelet of 6.5 /im thickness. 



amplitude. For spin that are at Co, the signal is propor- 
tional to cos ip. The data are fitted by the functional form 
exp(— Tp/r) cos(7£?iTp) + b. The period gives a calibration 
of the r.f. field strength at the sample location and we get 
B\ = 15G during the pulse. The other fitting parameters 
are r = 43 ± 6^.s and a positive offset b — 3 ± 0.2A. The 
values of these last two parameters depend strongly on 
B\ . The positive offset b is mainly due to the non-uniform 
field inside the sample [33]. For 1 H away from Co, the di- 
rection of B e is not exactly perpendicular to k and only 
a partial inversion of the z component is obtained when 
<p = tt. The decay of the magnetic moment fitted by r is 
due to field inhomogeneity which causes a dephasing of 
the magnetization in the transverse plane [34] . 

To study the transverse magnetization decay of 1 H [17] 
a sequence of 3 pulses is used. A n/2 pulse is applied to 
the spin system, so that the magnetization at Co is rotated 
to the transverse plane. After a fixed delay r , a tt pulse 
is applied to form a spin echo. Shortly after, a tt/2 pulse 
takes an instant snap-shot of the transverse magnetiza- 
tion by rotating it along k and the frozen component is 
measured with the c.w. sequence. Varying the time delay 
Tb between the last two pulses reconstructs the transient 
shape of the spin echo. Using the same settings as the ear- 
lier measurement, the widths of the 7r/2 and n pulse are 
set to 3.8/xs and 7.6/zs respectively. The delay between the 
center of the first two pulses is r a — 17/xs. In Fig. 5, the 
lock-in peak (again averaged over Is around its maximum) 
is shown as a function of r a + t&. As expected for a spin 
echo, the reconstructed transverse magnetization becomes 
refocused at a time 2r a . 




20 40 50 60 70 80 
time t. + r b (its) 



Fig. 6. Spin-spin relaxation time measurement: Normalized 
heights of the spin echo are displayed on a square-logarithmic 
scale as a function of r a + Tb with r a = Tb ■ The straight line 
is a fit with exp {-(2r a /T 2 ) 2 } where T 2 = 39 ± 1/js. 



With increasing spacing r between pulses, the size 
of the spin echo signal decreases due to spin-spin relax- 
ation. Using the same sequence as above, Fig. 6 is a plot 
of the spin echo amplitude measured as a function of the 
time 2r a . The measured values are plotted on a x 2 - log(y) 
scale and one finds that the data follow the relationship 
exp{-(t/T 2 ) 2 } with T 2 = 39 ± lfts. With the inferred 
T 2 , the shape of the echo in Fig. 5 can be calculated tak- 
ing into account the dipolar linewidth of the protons in 
our compound [28] and the spatial dependence ai(C). The 
solid line in Fig. 5 is the best fit obtained for a sample 
thickness of 6.5^im which is in good agreement with the 
value obtained on the image. 

The longitudinal magnetization recovery is now mea- 
sured after a saturation comb [35] . This protocol puts effi- 
ciently inhomogeneous spin systems in a well defined uni- 
form state outside thermal equilibrium. The saturation 
comb is composed of three n/2 pulses spaced by 100/xs. 
The c.w. sequence is applied at a variable delay (13ms 
< t < 20s) after the comb. In order to obtain an intrinsic 
measurement of the relaxation, it is important to ensure 
that the sensitivity profile a\ (C) is exclusively included in- 
side the sample section, otherwise a partial re-polarization 
of the magnetization occurs during the measurement cy- 
cle [17]. For our settings, Co is set exactly at the middle 
of the sample and r — 2.4/xm is smaller than the sam- 
ple thickness. As before, the value plotted is the lock-in 
output averaged over a Is time interval around its maxi- 
mum. No signals are detected when t = 13ms. On Fig. 7, 
two relaxation times in the recovery process are clearly 
observed. The results are fitted with a double exponential 
g s {1 - exp(-t/T ls )} + (1 - g s ) {1 - exp(-t/Ti,)} which 
gives g s = 49±2%, T ls = 0.35±0.03s and T u = 5.4±0.5s. 
The value g s does not correspond directly to the propor- 
tion of spin that have a short relaxation (n s ) since the 
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Fig. 7. Measurement of the longitudinal magnetization recov- 
ery: The logarithmic of the normalized amplitude of the force 
signal is shown as a function of the interval t between a satu- 
ration comb and the c.w. sequence. The solid line is a fit with 
a double exponential recovery which yields Ti s = 0.35 ± 0.03s 
and Tu = 5.4 ± 0.5s. Each point is the average of 32 c.w. se- 
quences. A r.f. field of Bi — 15G is used for the comb. The 
settings for the cyclic adiabatic inversions are a r.f. field of 7G, 
O/2-k — 50kHz and r; = 100ms. We make sure that there is no 
net repolarization of M during the c.w. sequence. 



factor between M z and the lock-in amplitude is also a 
function of T\ . 

Neutron diffraction studies [36] of the (NH^SC^ crys- 
tal structure show that there are two NH4 sites at room 
temperature surrounded respectively by five and six SO4 - 
ions. The protons of the two inequivalent ammonium ions 
are coupled via dipole-dipole interactions and the mea- 
sured spin-lattice relaxation rate at 300K is an averaged 
value of the . In a variable temperature NMR measure- 
ments, O'Reilly and Tsang [23] observe a single exponen- 
tial 1 H relaxation process and analyze their T\ results by 
the reorientation correlation times To of the two distinctive 
NH^ . At 300K, 1/to should be larger than the Larmor fre- 
quency and the rotation should be isotropic which means 
that Ti should be independent of the orientation between 
the static field and the crystallographic axis of our sam- 
ple. We suppose that our observed two T\ processes might 
be due to water contamination inside the sample during 
its contact with air. The presence of H 2 in the crystal 
lattice could decrease the reorientation correlation time 
of the ammonium ions, hence diminishing the protons T\. 
The relatively high proportion of spin with short 7\ might 
be due to the exceptionally small thickness of our crystal 
(7 fan). 

In order to check the later hypothesis, a conventional 
NMR measurement was performed by A. Dooglav with a 
IT custom spectrometer. The sample consisted of ~ lg of 
our sample ground to small particles with dimensions be- 
low 50/im. In this fine powder sample, a double relaxation 
process is also observed with the following parameters 
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Fig. 8. Force signal as a function of B ex f. A saturation comb 
is applied 0.6s (closed circles) and 16s (open circles) before 
the c.w. sequence. The solid line is the expected profile for a 
parallelepiped sample of 7/im thickness within both the free 
spin and adiabatic approximations. The c.w. sequence uses a 
r.f. field of 7G, H/2n = 25kHz and n = 300ms. The inset 
shows the transfer function that corresponds to these settings, 
01 is the spatial dependence of the sensitivity profile. 



n s = 17±5%, T ls = 0.37 ± 0.1s and T u = 4.7 ± 0.2s. This 
result is in sharp contrast with experiments performed on 
coarser grains, where only one relaxation is observed with 
T\ = 5.0 ± 0.2s. The values of the two relaxation times are 
equal, within error bars, to the ones measured by MRFM. 
In addition, measurements were performed on the same 
50^m powder after two weeks of aging in air. It showed 
a rise of n s to 26 ± 3% in the longitudinal magnetiza- 
tion recovery experiment. A standard spin-spin relaxation 
measurement on this powder seems also to indicate a dou- 
ble time T 2 with n s = 20 ± 6%, T 2s = 49 ± 12/xs and 
T 2 i = 79 ± 1/xs. 

One corollary issue concerns the spatial distribution 
of each spin species inside the sample section. To perform 
this measurement, we record the amplitude of the lock-in 
signal as a function of B ext for two delays t between the 
saturating comb and the c.w. sequence. By sweeping B ext , 
the surface Co is displaced to a different height in the sam- 
ple. The force signal is then proportional to the density of 
spin around this location. In order to obtain a local mea- 
surement, the thickness r of the slice probed is reduced by 
decreasing both fi and B\ during the c.w. sequence. The 
inset of Fig. 8 shows the spatial dependence of the transfer 
function ai(C) for our settings where r, the half width 
at half maximum, is 1.9/xm. By changing the delay t, the 
weight q s of one spin species compared to the other can 
be varied. Qualitatively, the measurement protocol gives 
more weight to the spin species with short relaxation when 
the comb is close to the c.w. sequence. The obtained re- 
sults are shown in Fig. 8 for both t = 0.6s (closed circles) 
and t = 16s (open circles). A first look at the result in- 
dicates that a more rounded distribution is obtained for 
the t = 0.6s data. The measurements, however, collected 
close to the edge of the sample are skewed by repolariza- 
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tion processes that modify the shape of the lock-in signal. 
Inside the bulk of the crystal (0.942T < B ext < 0.944T), 
there is no clear evidence of a spatial modulation of one 
spin population compared to the other, e.g. a dip of the 
signal in the middle of the crystal. This result suggests 
that, within our resolution, the water contamination is 
uniform in the thickness. The solid line is a calculation of 
the expected profile for a parallelepiped sample of dimen- 
sions 100 x 50 x 7/im 3 within both the free spin and adia- 
batic approximations. In spite of the idealized model, the 
t = 16s data (open circles) are well described by the calcu- 
lated profile except for the high field range. The shoulder 
at B ext — 0.946T corresponds to the surface of the sample 
that has been glued with epoxy to the cantilever. We did 
not attempt to fit this part of the data. The observed step 
in the signal might be due to the protons in the epoxy. 
A small roll angle between the sample and the cantilever 
combined with the particular shape of our crystal is also 
consistent with the observed effect. 

Although the data analyzed here-above ensure that 
two populations of spin with different NMR properties 
are present, their actual proportion is not quantitatively 
determined, as the actual values of the relaxation times in- 
fluence the magnitude of the measured signal. For a bet- 
ter knowledge of the sensitivity of the technique to the 
measurement parameters, it is then necessary to perform 
quantitative analyses. 



4 Quantitative measurements 

In this section, we shall first calibrate the mechanical re- 
sponse of the cantilever and the mechanical noise. The 
time dependence A(t) of the lock-in signal is calculated, 
taking into account relaxation processes and non-adiabatic 
effects. The experimental responses for different values of 
Bi and Q are compared with the calculations. This al- 
lows us to select an experimental condition for which non- 
adiabatic effects can be neglected. It is then shown that 
two relaxation times are indeed required to fit the time 
dependence of the observed lock-in signals, with values 
consistent with those obtained from T\ data. 

The Q of the cantilever is first measured carefully throu- 
gh the noise vibration spectrum of the cantilever loaded 
with the sample in vacuum when the r.f. power is off. The 
lock-in time constant is set to 77= 10s. An audio generator 
sweeps the lock-in reference around u c . The plotted value 
in Fig. 9 is the standard deviation of the lock-in signal 
estimated over 100 xt; (the mean lock- in signal is zero). 
During the whole experiment, the temperature stability 
of the cantilever is better than ±0.01°C which guarantees 
that uj c does not shift by more than 0.01Hz. Fitting the 
squared amplitude with a Lorentzian [37] , one obtains the 
cantilever resonance frequency lu c /2tt = 1397.77Hz and 
quality factor Q = 4000 (defined as the ratio of uo c over 
the full width at half maximum of the power spectrum). 
Away from resonance, our sensitivity is limited by the 
noise of the detection electronics. It is several orders of 
magnitude smaller than the A-scale motion of the can- 
tilever at resonance and therefore it can be neglected. 
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1399 1400 



Fig. 9. Noise vibration spectrum of the cantilever with the 
sample attached: The data is obtained in vacuum (10" 2 torr) 
when no r.f. fields are applied. The lock-in time constant is r; = 
10s. The signal is the standard deviation of the lock-in output 
estimated over 100 xti. The cantilever holder temperature sta- 
bility is better than ±0.01° C during the measurement. 



Near lo c , the cantilever motion consists of white noise am- 
plified by a narrow-bandwidth mechanical resonator [38]. 
Av c is the one-sided equivalent noise bandwidth (ENBW) 
of the mechanical resonator Av c — uj c /8Q — 0.27Hz. The 
noise at the output of the lock-in is this narrow-band mo- 
tion noise observed through a RC filter of time constant 
n = l/RC whose ENBW Av t = 1/4-rj = 0.025Hz. Exactly 
at resonance, the combined distribution gives an ENBW 
Av = {l/Av c + l/zA^) -1 =0.023 Hz. To convert our data 
to spectral densi ty, t he resonance amplitude in Fig. 9 must 
be divided by V 'Av. The measured noise spectral density 
is A = 13A/\/Hz. This figure also corresponds to the noise 
observed in Fig. 3 where Av = 0.20Hz. The result has to be 
compared with the intrinsic correlation function for fluc- 
tuations of a Brownian particle harmonically bound to a 
single degree of freedom oscillator of spring constant k: 



(AA(t)) = 



4Qk B T 

kljJr 



(5) 



Taking the square root of the above expression, one gets 
A T = lOA/VHz at T = 300K which is in good agree- 
ment with our measured value. In conclusion, our domi- 
nant noise comes from the thermal vibration of the can- 
tilever. From this result, one can estimate the smallest 
force detectable by the instrument in one shot kA/Q = 

2xlO- 15 N/ v / Hi. 

In order to obtain a quantitative measurement of our 
force signal when the r.f. field is applied, a more detailed 
study of the time dependence of the lock-in signal A(t) 
is needed. The length of the c.w. sequence is increased to 
6s compared to Fig. 3. Fig. 10 is the average of the lock-in 
signal over 32 sequences using a short lock-in time con- 
stant of T[ = 30ms. The striking feature of this plot is 
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time (ls/div) 

Fig. 10. Time dependence of the lock-in signal: The trace is 
the average of 32 c.w. sequences with t; = 30ms. The c.w. se- 
quence is composed of a r.f. field of 10G with i?/27r = 50kHz. 
The solid line is the calculated vibration amplitude of a har- 
monic oscillator driven by a damped sinusoidal magnetic force 
with a decay time constant r m = 2.2 ± 0.07s. In this case, the 
predicted maximum signal is ^peak = 0.5Aq with Ao = QFo/k. 



ation mechanisms in this compound are associated with 
the time varying field induced by the change in the NH^ 
orientation. For an exponential correlation function with 
correlation time To, T\ x is expressed as a sum of the spec- 
tral density of these fluctuating fields J^(lu) = t /(1 + 
co 2 tq) with an index i that corresponds to the number of 
net spin flip: 1/T lx = 3/2 7 4 ft 2 /(7 + l)/r 6 (5/2j( 1 )(w ) + 
1/4J( 2 )(2w ) + l/4j(°)(2w e )), with uj e = lB e . In the ap- 
proximation for our compound that w e To <C 1, one ob- 
tains that Ti x is independent of 9 and T\ y = T\ x . Coming 
back to equation (6), one recognizes that the first inte- 
gral is due to the finite value of the non-adiabaticity pa- 
rameter a (see appendix B). The second represents the 
spin-lattice relaxation T\ p in the rotating frame [41]. The 
last integral is the equilibrium magnetization that corre- 
sponds to the spin temperature in the rotating frame, r 
and M eq are defined by rewriting the above expression in 
the form Z\/27r = —M/ (uj c t) + M cq / '(lj c t) . During the c.w. 
sequence, the oscillatory driving magnetic force is damp- 
ened and the instantaneous value is given by: 

F(t)mg sin(uj m t) / a x (C) M((, t) d( + constant, (7) 



that the norm of the magnetization | M | decays dur- 
ing the c.w. sequence. We have checked that experimental 
perturbations such as the phase noise of the r.f. source are 
negligible in this case [39]. 

In the rotating frame, the magnetization tends to re- 
cover slowly towards a steady state value due to spin- 
lattice relaxation processes [40] . These relaxation processes 
are different from the relaxations measured in section 3 
which are linked to the time dependence of the magneti- 
zation in the absence of r.f. fields. In addition, the lack 
of adiabaticity quantified by a w 9/"/B e produces a pre- 
cession movement around the locally changing effective 
field direction. This mistracking corresponds to a mag- 
netization component perpendicular to the instantaneous 
precession axis B e + 9/jj that relaxes due to spin-spin in- 
teractions. But in the limit where a < 1 and 2-7T /uj c <C T\, 
the decrease of M after one cycle A = M(t+2ir/u) c )—M(t) 
is small. The value is calculated at the lowest order for one 
spin species (see appendix B). 



with 




x cos ( [ 7 f,,^ dt''\ dt'dt 
\J t , sin 0(f) J 

cos 2 8(t) sin 2 9(ty 



+M 



s: 



cos 9(t) 



lz 



dt, 



'- lx 



dt 



with Ti z is the usual Xi in the absence of r.f. field, T\ x is 
the transversal spin-lattice relaxation and 1/T X + = (1/Ti y + 
cos 2 9/Ti x + sin 2 9/Ti z )/2. It was shown that the relax- 



M((,t) = \ M eq (C) + {M - M oq (C)}exp 



-(C) 



(8) 



The integral in equation (7) relaxes approximately accord- 
ing to a single exponential towards its equilibrium value 
rn cq = f v ai(()M cq (()d( with an apparent characteristic 
time r m . One notes that m oq = by symmetry when £o is 
centered at the middle of the sample (o = £ TO . The value of 
TO C q is positive for Co < (m and changes sign for Co > Cm 
[17]. In the particular case where m oq = and a <C 1, 
then it can be shown that the coefficient r m is bounded 
between T lx < r m < T lz [21]. 

The forced vibrations of an harmonic oscillator are 
given by the convolution product: 



a(t) =/3 [ 
Jo 



4 F(t') f t-t' 
— cxpi 



sin {to c {t — t')} lo c dt', 
(9) 



with (3 = {l + 1/(4Q 2 )} and l/r c the damping constant 
of the cantilever. In our experiment the external force 
F(t) — Fq exp(— t/T m ) exp(iuj m t) for m oq = 0, with 



is 



Fq = kAo/Q and r TO the characteristic decay time of the 
magnetic force. a(s) the Laplace transform of equation (9) 
is calculated in the complex plane: 



k 



(6) TT«(s) 



Tm (4Q 2 + 1) 



(sr m + 1 - ir m u m ) {(st c + l) 2 + 4Q 2 } ' 



(10) 



The response A{t) of the lock-in is the imaginary part of 
the inverse Laplace transform {a(s + iuj m )/(l + r;s)} 
with t i the lock- in time constant. An approximation can 



O. Klein, V. V. Naletov, and H. Alloul: Mechanical detection of nuclear spin relaxation in a micron-size crystal 



9 





50 r 




40 




30 - 




3 

"S. 


20 


s 


10 - 








(a) experiment 

/ t 



„° o 14G 



X 



10G 



-0 7G 



50 100 150 200 250 300 
n/2n (kHz) 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




1 1 1 

(b) theory 


i i i 




^ - 14G 








10G 




^ 7G 







50 100 150 200 250 300 
Q/Z-n (kHz) 



Fig. 11. (a) Measurements of the amplitude of the lock-in sig- 
nal is shown as a function of the modulation amplitude for dif- 
ferent strength of the r.f. field. The dashed lines are guides for 
the eye. The arrows indicate the limit of the adiabatic regime. 
The measurements are performed at B ex t = 0.9425T. (b) Cal- 
culated amplitude of the lock-in signal obtained for a 7/im thick 
sample using equation (7). The parameters introduced in the 
model are T lz = 4.9s, T lx = T ly = 3.2s and T 2 = 40^s (solid 
line) or 100/is (dotted line). 



be obtained in the special case where cu m — uj c in the limit 
Q > 1 and uj c > (l/r m , 1/n). 



A a 
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-n) 



exp l - — 



exp 

Tl 



exp [ - — 

' m. 



(11) 



The sum of these three exponentials vanishes at t = and 
each term decays to zero with a different time constant 
at later time t > 0. This leads to a peaked lock- in signal 
^peak whose amplitude and position depends on r m (for a 
fixed r c and n). At the end of the c.w. sequence in Fig. 10 
(t > 6s) , the free oscillations decay of the cantilever (time 
constant r c ) are observed. If one tries to fit the data with 
the above nonlinear form, r m = 2.2 ± 0.07s is obtained but 
the quality of the fit is not very good. Values of r m smaller 
than Ti x = 3.2s have also been reported by Verhagen et 
al. [21] and these findings were attributed to the phase 
noise of the r.f. source. However, when large modulation 
amplitude are employed for the c.w. sequence, such a fast 
force decay can also be consistent with a magnetization 
decrease due to a lack of adiabaticity. 

To understand further the meaning of this fit param- 
eter r m , we plot in Fig. 11a the lock-in peak amplitude 



measured for different values of f2 and B\ when B ext — 
0.9425T. The value of the non-adiabaticity parameter a 
increases along the abscissa axis. For a fixed Q and B\, the 
value of a oscillates with time and passes through a maxi- 
mum, a max = fiuJm/^Bl, at time a t = modulo ir /uj m . 
Fig. lib shows the amplitude of the peak signal A pca ^ pre- 
dicted by equation (11) with F(t) calculated from equa- 
tion (7) using a sample thickness of 7/im. The results are 
normalized by A tot = QgMoV s /(ky/2) the amplitude as- 
sociated with a uniform inversion of all spin inside the 
sample. The parameters introduced in the model are the 
values of the spin-lattice relaxation times measured on 
powder samples by conventional NMR with T\ z = 4.9s 
along the static field and T\ x = 3.2s along a 10G r.f. field, 
in the approximation that T\ y = T\ x . In our theoretical 
model, perturbation effects from the dipolar broadening 
are also introduced. They are approximated as a time in- 
dependent local field of Lorentzian lineshape. The solid 
lines are the shape calculated with T 2 = 40/^is (the dotted 
lines correspond to T 2 = 100^<s). 

The increase of the force signal at small ]?(< 50kHz) 
corresponds to an increase of the modulated magnetic mo- 
ment. A larger frequency deviation increases the width of 
the probed slice, -T, and more protons oscillate at the fre- 
quency uj m . For both large B\ > I4G and large Q > 
1 50kHz the amplitude of the signal eventually saturates 
when r becomes larger than the sample thickness. 

For B\ = 7 and 10G, the deviation from the low fl- 
lincar increase (indicated by the arrows) marks the cross- 
over from an adiabatic regime to a quasi-adiabatic one 
[40]. In our sample the threshold occurs at a max = 0.1, in 
good agreement with the theoretical model. In the adi- 
abatic regime (a max < 0.1) the value of r m is deter- 
mined by Ti p effects, while in the quasi-adiabatic regime 
(0.1 < a max < 1) the magnetization decay r m becomes 
smaller than T\ x . The predicted position of this cross-over 
depends somewhat on the shape of the proton lincwidth 
(dotted lines). 

From the last discussion, one can conclude that the 
settings of Fig. 10 correspond to a non- adiabatic parame- 
ter a max = 0.04 well inside the adiabatic regime for our 
compound. It can then be inferred that the decrease of 
the force signal in Fig. 10 is due to spin-lattice effects in 
the rotating frame and our fitted value of r m must be an 
average of the two T\ reported in Fig. 7. 

The time dependence of the lock-in signal is fitted 
with a double damped synchronous excitation of two spin 
populations with respectively short (r ms ) and long (r m i) 
relaxation times. The nonlinear function n s A(t, r ms ) + 
(1 — n s )A(t, T m i) is used where A{t) is given by equa- 
tion (11). The best fit is obtained for r ms = 0.55 ± 0.02s, 
T m l — 4.8 ± 0.2s and n s = 69 ± 1%. The result is the solid 
line shown in Fig. 12. The fit values obtained for the relax- 
ation times are similar to those measured in the magneti- 
zation recovery experiment in Fig. 7. On Fig. 12 the sepa- 
rate contribution to the lock-in signal of each spin species 
(dashed lines) are shown. The maximum force signal of 
the two spin species occurs respectively 0.7s and 1.9s af- 
ter the start of the c.w. sequence for the short and long 
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Fig. 12. Signal of Fig. 10: The solid line is the response of 
an harmonic oscillator driven by two synchronous forces with 
short and long relaxations times n s A(t, T ms ) + (l — n s )A(t, r m i). 
The best fit is obtained for r ms = 0.55 ± 0.02s, r m ; = 4.8 ± 0.2s 
and n 3 = 70%. The height of the peak is 0.3Ao. The dashed 
lines show the contribution of each spin species. 



T m . It brings up the question on how to define the lock- in 
peak amplitude in the case of a sample containing two spin 
species. We recall that our definition of the force signal is 
the average of the lock-in peak amplitude over a Is time in- 
terval around its maximum value. This approach gives ap- 
proximately equal weights to both spin species in the mea- 
surement. One can also observe in Fig. 12 that the height 
of the two peaks are approximately equal despite the fact 
that there is 2.3 times more spin with short relaxation. As 
a matter of fact it can be shown that the mechanical de- 
tection is 2.4 times more sensitive to spin that have a 4.8s 
relaxation time compared to spin that have a 0.55s one. 
From this result, the value of the fit parameter g s w 0.5 
in Fig. 7 can be converted into the proportion of spin that 
have a short relaxation n s — 2.4g s /{(2.4 — l)g s + 1} = 70 
%; a value that agrees well with the fit n s in Fig. 12. 

Finally, the expected amplitude of the force signal for 
our sample is calculated. Using equation (4), one gets a 
value of A$ = 100A for the settings of the c.w. sequence 
used in Fig.12 (B 1 = 10G and S7/2ir = 50kHz). In Fig.12 
the predicted maximum of the lock-in peak is 0.3.Ao = 30A 
which is close to the experimentally measured value of 
20A. In conclusion, our measured amplitude of the peak 
lock-in signal is in good agreement with the theoretical 
prediction if the two spin-lattice relaxation times of the 
two spin species are taken into account. Other effects such 
as misalignment of the sample compared to the cylinder 
magnetic axis can account partially for a decrease of the 
signal (e.g. a small offset of 0.1mm from the axis decreases 
the amplitude of the lock-in signal by a factor of 2). 



5 Conclusion 

Measurement sequences combining fast adiabatic passages 
and pulses are reported. They allow us to measure T\ and 



T2 for microscopic samples using a mechanical detection. 
This has been applied to quantitative analyses of the de- 
tected signals for a 7/^m thick sample of (NH 4 ) 2 S04. The 
transverse relaxation T2 has been found consistent with 
conventional NMR detection on a macroscopic sample. 
Our sample displays however two spin lattice relaxation 
times Tis = 0.4s and Ty = 5s. While the long Ty corre- 
sponds to that measured for coarse powder samples, the 
short Ti s might be due to water contamination of our 
7/xm thick crystal during its contact with air. This con- 
tamination is found to be uniform in the thickness of the 
sample. This large difference in T\ values has allowed us 
to study the influence of the spin-lattice relaxation in the 
rotating frame on the measured time dependence of the 
lock-in signal, as well as the variation of signal intensity 
with increasing non-adiabaticity of the sweep sequence. A 
consistent analysis of all the experimental parameters has 
been proposed and will be quite useful in future quan- 
titative investigations of MRFM signals. Our work also 
raises the problem on how to perform reliable spin lat- 
tice relaxation measurements at different locations in the 
sample. Our investigation is mainly restrained to the bulk, 
i.e. the middle of the sample. It is known that the time 
dependence of the lock-in signal (and thus the apparent 
relaxation times) depends strongly on the value of B ext . 
These difficulties prevented us from interpreting quantita- 
tively our results on the spatial distribution of the different 
spin densities close to the sample surface. This issue will 
be best resolved by performing a similar experiment on a 
hetero-layer sample of well-characterized composition. 
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A Inhomogeneous field 

Near the axis, a uniformly magnetized (M a ) cylinder of 
length I and diameter <j> produces a field, whose component 
along fc, B cy i 7 decays radially as 



B cy i{r,z) 



4ttM, 



61 



+3 6 



z + l 



z + 1 



bi 



+ G(r 4 ), (A.l) 



with b a (z) = z (l + 4z 2 ) ° . The fields are expressed in 
cylindrical coordinates with the origin centered on the 
cylinder upper surface (see Fig.l). In our case M s is cal- 
culated from the applied field B ext needed to produce a 
resonance signal at the sample position, z = 0.70mm. 
B cy i(0,z = 0.70) = wo/7 - B ext = 0.352T is replaced in 
the above expression and one obtains M s w 1400emu/cm 3 
for our iron. Using this result, the gradient g = — 470T/m 
at the sample location can be calculated. 
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B Adiabaticity 

The aim of this appendix B is to calculate the decrease of 
the magnetization due to the spin-lattice relaxation and 
the lack of adiabaticity. The solution below is proposed by 
M. Goldman. In the limit of strong r.f. fields (larger than 
the local field), one can neglect the spin-lattice relaxation 
of the dipolar energy expectation value. In the rotating 
frame, the time evolution of the different spin components 
are [40]: 



d{h) „ , M , (Io)-(h) 



dt 



(B.la) 



dt 
d(I ; 



+ 7 B 1 cot 8(I y ) - - U[U Dz ,I x ]) (B.lb) 
J-ix n 

'" = + 7 B 1 {{!„) - cot 6(I X )) - £!L - l([H Dz m$q) 



dt 



Tly k ' 



with (I) = Tr(Jer) the expectation value of the magneti- 
zation, er the instantaneous density matrix in the rotating 
frame and TLdz the secular part of the dipolar Hamilto- 
nian. The commutator incorporates the local field contri- 
bution defined through B\ = Tr(H 2 Dz ) / 7 2 Tr(7f). In our 
notation 8 is the angle between the directions of the static 
and effective field, B e = B\j sin#, with B\ cot# the pro- 
jection along k. Under r.f. irradiation, a new coordinate 
system {A, Y, Z} is defined through a transformation by 
the unitary operator exp(—i0I y ), a rotation around y by 
9. In the doubly rotating frame the differential equations 
then become: 



s 2 c 2 



d(Iz) (Io) , • , 

~dT~ i" ' 7^ 

io-csU--^}: / 



(B.2a) 



d(Ix) 
dt 



1 1 

T\ x T\ z 



(Iz) 



B\ . . f c 2 s 2 . , , 
+7 — (I Y )- i — + — }(Ix) 

S I -tlx J-lz 



3c 2 -l i 
2 h 



[HdztIx]) 



(B.2b) 



- -l—^X) - — 7, T{[rLDZ,lW)t c ) 



dt 



-1» 



2 h 



with s = sm9 and c = cos6. Hdz is the doubly truncated 
part of the dipolar Hamiltonian that commutes with Iz- 
The term 6/jj is the inertia field due to the transforma- 
tion to a time-dependent reference axis. In the adiabatic 
regime, defined by a = 8/(-f B e ) <C 1, it can be shown that 
(Ix) = (Iy) = and the first two terms of equation (B.2a) 
are the expression of the spin-lattice relaxation in the ro- 
tating frame for strong r.f. fields [1,41]. This appendix 
seeks to evaluate the term proportional to (Ix) in equation 
(B.2a) that represents the decrease of (Iz) due to the lack 
of adiabaticity. Our investigation will be restricted to the 
quasi-adiabatic regime [40], where non-adiabaticity cor- 
rections might be dominant over T\ p effects but a < 1. The 



equations of motion are expressed in terms of the raising 
and lowering operators (I + ) = (Ix) + i(W) and (I~) its 
complex conjugate. We suppose that -fB e <C 1/t , the re- 
orientation correlation rate, for our compound which leads 
to T\ x independent of 9 and T\ y ss T\ x . As a further sim- 
plification, spin-spin interactions are neglected and these 
effects will be calculated numerically later on. As a con- 
sequence, one has 



d(I + ) 
dt 



- o 



cs 



\T lx 



1 

T\ z 



{Iz) 



J-u 



with 1/T+ = (1/T ly + c 2 /T lx + s 2 /T lz )/2. Furthermore, 
the period of the cyclic passage is much smaller than the 
spin-lattice relaxation times. Hence, both 1/T\ Z and 1/Tia, 
are negligible compared to 9: 

(B.4) 

which, upon integration, gives the result: 



(I + ) = - £ ! [ 9(t')(Iz)+s\n9(t')^L 



x cxp 



assuming that (I + ) = at a time t = 0. The expression 
for (I~) is the complex conjugate of the above expression. 
For values of t < Ti, it is a good approximation to ne- 
glect (Iq)/Ti z compared to 0(lz) in the first bracket. Since 
the decay of (Iz) is slow, it is replaced by a constant. Fi- 
nally, the time variation of the longitudinal magnetization 
is given by: 



d(Iz) 
dt 



x cos 



7*1 - 



-(Iz) 
+ (Io) 



t , sin6»(t" 
cos 2 0(t) sin 2 9(t) 



dt" dt' 



T u 
cos9(t) 
~Tu~' 



-lx 



(B.6) 



which has been used in equation (6) in the text. Although 
the general form of our final expression contains a T+- 
exponential decay in the first integral, the decrease of 
the magnetization due to the quasi-adiabatic part of the 
c.w. sequence is unrelated to spin-lattice relaxation mech- 
anisms. This is best seen in the free spin limit (Ti — ► oo), 
where the first term in equation (B.6) reduces to the inte- 
gral form of the ordinary differential equation d(I) /dt = 
j (I) x B e , an equation of motion that preserves the norm 
of the magnetization. Our above solution expresses the 
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decrease of (Iz) with the set of initial condition (I + ) = 
at t = 0. Without local field effects, the free spin solu- 
tion would eventually reach a finite "steady state" and for 
these large time scales, it is then important to take into 
account spin-spin interactions. 

These interactions are best seen in another rotating 
frame {X 1 , Y' , Z'} where the direction Z' is aligned along 
the instantaneous axis of precession of the magnetization. 
The new system of differential equations is obtained by 
applying the unitary transformation exp(— ialx) to the 
equation system (B.2). Here second order corrections are 
examined and a more detailed analysis should take care of 
the new inertia field a. In this triply rotating coordinate 
system, the relaxation mechanisms along the X' and Y 1 di- 
rections are dominated by local field effects and these have 
a characteristic time of the order of Ti ~ 40^.s in our sam- 
ple. The exact expression is more complicated because the 
local field in the triply rotating frame oscillates with time. 
The system of differential equations is solved numerically 
in the approximation that the local field is time indepen- 
dent and of Lorentzian lineshape. The result is compared 
with the expression obtained in equation (6). The two cal- 
culated values of (Iz 1 ) (~ (Iz) when a < 1) are close in 
both cases after one passage (n/w m ). The main difference 
occurs after several passages: when local field effects are 
included, the amplitude (Iz 1 ) decays toward zero and the 
expectation value of the norm of the magnetization follows 
approximately the instantaneous value of (Iz 1 ) ■ 
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